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MINIMUM-MASS  DESIGN  OF  FILAMENTARY  COMPOSITE  PANELS 
UNDER  COMBINED  LOADS:  DESIGN  PROCEDURE  BASED 
ON  SIMPLIFIED  BUCKLING  EQUATIONS 

W.  Jefferson  Stroud  and  Nancy  Agranoff 
Langley  Research  Center 

SUMMARY 

An  analytical  procedure  is  presented  for  designing  hat- stiffened  and  corrugated 
panels  made  of  composite  material  and  subjected  to  longitudinal  (in  the  direction  of  the 
stiffeners)  compression  and  shear  loadings.  The  procedure  is  based  on  nonlinear  mathe¬ 
matical  programing  techniques  and  a  simplified  set  of  buckling  equations.  Design  require¬ 
ments  considered  are  buckling,  strength,  and  extensional  and  shear  stiffness.  Studies 
show  the  effects  of  specified  thickness,  variation  of  cross-section  dimensions,  stiffness 
requirements,  local-buckling  boundary  conditions,  and  the  effect  of  combined  compression 
and  shear  loadings. 


INTRODUCTION 

Because  composite  materials  offer  the  potential  for  substantial  mass  and  cost 
savings,  Langley  Research  Center  is  engaged  in  a  broad  program  to  study  structural 
applications  of  filamentary  composite  materials.  One  objective  of  that  program  is  to 
establish  a  mass  and  strength  data  base  for  generic  structural  components  made  of  com¬ 
posite  material.  Such  a  data  base,  containing  both  analytical  and  experimental  results,  is 
necessary  to  develop  confidence  in  composite  structural  design  procedures.  The  generic 
component  considered  in  this  report  is  the  stiffened  panel:  specifically,  hat-stiffened  and 
corrugated  panels. 

Some  of  the  previous  work  dealing  with  the  analysis  and  design  of  stiffened  composite 
panels  is  discussed  in  references  1  to  6.  The  present  report  presents  an  analytical  pro¬ 
cedure  for  designing  low-mass,  hat-stiffened  and  corrugated  panels  which  are  subjected  to 
longitudinal  compression  and  shear  loadings.  Design  studies  carried  out  with  the  pro¬ 
cedure  are  also  presented. 

The  procedure  is  based  on  nonlinear  mathematical  programing  techniques  in  which 
buckling,  strength,  and  stiffness  requirements  as  well  as  upper  and  lower  bounds  on  cross- 
section  dimensions  are  taken  into  account  with  inequality  constraints.  The  objective  is  to 


obtain  minimum-mass  designs.  Overall-buckling  loads  are  calculated  by  assuming 
smeared  stiffnesses.  Local-buckling  loads  are  calculated  by  assuming  simple -support 
boundary  conditions  along  the  lines  of  attachment  of  adjacent  elements.  Certain  buckling 
modes  (stiffener  rolling,  for  example)  are  neglected  in  the  buckling  analyses  used  herein. 
Lamina  stresses  and  strains  are  used  in  the  strength  calculations. 

In  composite-panel  designs,  the  interaction  of  design  variables  may  be  quite  com¬ 
plex.  To  demonstrate  typical  design  study  options  available  with  the  procedure  presented 
in  this  report,  studies  were  performed  on  hat-stiffened  and  corrugated  graphite- epoxy 
panels  loaded  in  compression  and  shear.  These  studies  show  the  effects  of  discrete 
thickness,  variations  of  cross-section  dimensions,  stiffness  requirements,  local-buckling 
boundary  conditions,  and  the  effect  of  combined  compression  and  shear  loadings.  Detailed 
study  results  are  presented  so  that  a  designer  may  obtain  a  better  understanding  of  the 
mass  penalties  associated  with  practical  design  constraints  such  as  stiffener  cap  widths 
and  stiffener  spacing. 


SYMBOLS 

Values  are  given  in  both  SI  and  U.S.  Customary  Units.  The  calculations  were  made 
in  U.S.  Customary  Units. 


surface  area  of  panel,  pL 


area  enclosed  by  hat  (See  eq.  (A25).) 


laminate  stiffnesses  defined  in  equations  (Al)  and  (A2) 


value  of  for  element  k 


a 

^k 

D 


D. 


1] 


k 


panel  width  (See  fig.  5.) 


element  lengths  defined  in  figure  2 


depth  of  element  2  defined  in  figure  23 


laminate  stiffnesses  defined  in  equations  (Al)  and  (A3) 


value  of  Djj  for  element  k 


2 


Di,D2,D3, 


smeared  orthotropic  stiffnesses  defined  in  equations  (AlO)  to  (A13) 


^x’^y’^xyj 


d 

E 

El 


ET 

ETj^ 

Ei,E2 


G 

GT 

^12 

L 

Nx 

^XjCrit 

Nx 

Nxy 

Nxy,crit 

Nxy, 

Ny 


distance  between  stiffeners;  2b^ 

Young's  modulus 

beam  bending  stiffness  of  one  period  of  panel  cross  section;  defined  in 
equation  (A19) 

longitudinal  extensional  stiffness  of  panel;  defined  in  equation  (A5) 

longitudinal  extensional  stiffness  for  element  k,  defined  in  equation  (A4) 

Young's  modulus  of  composite  material  in  fiber  direction  and  transverse  to 
fiber  direction,  respectively 

shear  modulus 

shear  stiffness  of  panel;  defined  in  equation  (A6) 

shear  stiffness  of  composite  material  in  coordinate  system  defined  by  fiber 
direction 

panel  length 

applied  longitudinal  compressive  loading  per  unit  width  of  panel  (See  fig.  1.) 

value  of  compressive  load  that  causes  buckling 

longitudinal  compressive  loading  in  element  k  (See  fig.  6.) 

applied  shear  loading  per  unit  width  of  panel  (See  fig.  1.) 

value  of  shear  load  that  causes  buckling 

shear  loading  in  element  k  (See  fig.  6.) 


applied  transverse  (y-direction)  loading  per  unit  width  of  panel 


Zjj  vertical  distance  from  neutral  axis  of  element  3  to  neutral  axis  of  element  k 

a  angle  element  2  makes  with  the  horizontal  (See  fig.  23.) 

shear  buckling  parameter 

^xy  average  overall  shear  strain  of  panel,  defined  by  equation  (A7) 

%’^y’^xy  strains  in  X-Y  coordinate  system 

,y  laminate  strains  for  element  k 
^k  ^k 

lamina  strains  in  coordinate  system  defined  by  fiber  direction 

1  A  12 

9  shear  buckling  parameter  defined  in  equations  (20)  and  (25),  and  fiber 

orientation  angle  shown  in  figure  1 

ji  Poisson's  ratio 

Poisson's  ratios  associated  with  bending;  defined  by  equation  (AlO) 
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^12’^21  Poisson's  ratios  of  composite  material  in  coordinate  system  defined  by  fiber 
direction 

p  density 

Q'l,a2,Ti2  lamina  stresses  in  coordinate  system  defined  by  fiber  direction 

T  shear  stress 

Superscript: 

a  allowable 

Subscript: 

k  integer  denoting  element  number 

LOADINGS  AND  CONHGURATIONS  CONSIDERED 

Loadings 

Two  loadings  are  considered:  a  uniform,  longitudinal  compressive  loading  Njj  in 
the  direction  of  the  stiffeners  and  a  uniform  shear  loading  N^y.  These  loadings  are 
shown  in  figure  1. 


Figure  1.-  Loadings  considered. 
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Configurations 


The  two  panel  cross-sectional  configurations  considered  in  this  report  are  shown 
in  figure  2.  The  cross-section  geometry  is  defined  in  terms  of  the  four  dimensions  bj, 
b2,  bg,  and  b^.  The  quantities  bj,  b2,  and  b3  are  the  widths  of  elements  1,  2,  and  3, 
respectively.  The  quantity  is  the  half -width  of  element  4.  Because  of  symmetry, 

*^3 

b^  is  replaced  by  in  the  corrugated  panel. 

The  panels  are  assumed  to  be  made  of  a  filamentary  composite  material  with  ply 
orientation  angle  9  measured  from  a  line  drawn  in  the  direction  of  the  stiffeners  (the 
x-direction)  as  shown  in  figure  1.  In  the  examples  presented  in  this  report,  only  two 
fiber  orientations  are  considered:  0=0°  and  6  =  ±45°.  The  45°  plys  are  assumed  to 
be  balanced  and  symmetric  and  are,  therefore,  denoted  ±45°. 


p 


1 


(a)  Hat- stiffened  panel. 

Figure  2.  -  Configurations  and  design  variables.  Element  identification 
numbers  are  indicated  in  figure  2(a). 
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—  p 


(b)  Corrugated  panel. 


and 


0=±45 


0=0° 


(c)  Filament  orientation  pattern. 
Figure  2.-  Concluded. 


All  elements  contain  ±45®  material.  Elements  3  and  4  can  also  contain  0®  material. 
The  configurations  are  defined  so  that  the  0®  material  and  the  ±45®  material  have  specific 
locations.  In  both  configurations  t3(0®)  refers  to  the  thickness  of  one  of  the  two  equal 
0®  layers  in  element  3;  t^(0®)  is  the  thickness  of  a  0®  layer  sandwiched  between  ±45® 
skins  in  element  4. 

The  hat-stiffened  configuration  is  a  uniformly  thick  corrugated  sheet  of  ±45®  material 
attached  to  an  equally  thick  flat  sheet  of  ±45®  material  with  0®  material  at  specified  loca¬ 
tions  in  elements  3  and  4.  The  corrugated  configuration  is  a  uniformly  thick  corrugated 
sheet  of  ±45®  material  with  0®  material  at  specified  locations  in  elem_ent  3.  In  both 
configurations,  the  thickness  of  a  ±45®  sheet  is  denoted  t2(±45®)  or  simply  t2*  The 
thicknesses  of  the  ±45®  sheets  are  taken  to  be  uniform  to  provide  a  design  that  is  relatively 
easy  to  fabricate. 


ANALYSIS-DESIGN  PROCEDURE 

Nonlinear  mathematical  programing  techniques  are  used  to  adjust  the  element  widths 
and  thicknesses  (design  variables)  to  provide  minimum-mass  designs  that  meet  the  design 
requirements.  The  design  requirements  (constraints)  considered  are  strength,  buckling, 
stiffness,  and  upper  and  lower  bounds  on  the  cross-section  dimensions.  The  objective  of 
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this  section  is  to  present  and  discuss  briefly  many  of  the  expressions  needed  to  evaluate 
the  strength  and  buckling  constraints.  Equations  for  the  extensional  and  shear  stiffnesses 
are  given  in  the  appendix  together  with  detailed  expressions  for  other  quantities  that  are 
used  in  this  report. 


Stress  Analysis 

The  longitudinal  and  transverse  strains  and  the  longitudinal  stress  resultants  are 
calculated  for  an  applied  longitudinal  compressive  load  Nx-  The  transverse  loading  Ny 
is  assumed  to  be  zero  for  the  overall  panel  and  for  each  element  that  makes  up  the  panel. 
The  shear-stress  resultants  and  shear  strains  are  then  calculated  for  an  applied  shear 
loading  Nxy-  Finally,  by  use  of  coordinate -transformation  and  stress-strain  relations, 
the  strains  and  stresses  are  calculated  in  the  directions  parallel  and  perpendicular  to  the 
fiber  direction  in  each  lamina  of  each  element.  All  elements  are  assumed  to  be  specially 
orthotropic  (ref.  7)  so  that  the  inplane  problem  in  each  element  is  given  by  equation  (A2) 
in  the  appendix. 

The  stress  analysis  equations  are  summarized  in  table  I  and  are  discussed  in  the 

text. 


Longitudinal  compression.-  All  elements  of  the  panel  cross  section  are  assumed  to 

have  the  same  longitudinal  strain  given  by  equation  (1).  An  expression  for  the 

longitudinal  extensional  stiffness  ET  of  the  panel  cross  section  is  given  in  the  appendix. 

Since  the  transverse  loading  is  assumed  to  be  zero,  the  transverse  strain  e  in  the  kth 

^k 

element  is  given  by  equation  (2).  The  longitudinal  loading  Nx  in  the  kth  element  is 

k 

proportional  to  the  longitudinal  extensional  stiffness  ETj^  of  that  element  and  is  given 
by  equation  (3). 


Shear. -  In  the  case  of  a  hat-stiffened  panel,  the  applied  shear-stress  resultant  N^y 
is  assumed  to  be  distributed  over  the  panel  cross  section  in  the  following  manner.  Ele¬ 
ment  4  carries  the  entire  applied  shear  stress.  In  the  region  of  the  hat  stiffener,  the 
applied  shear  stress  divides  into  two  parts.  Element  1  carries  one  part  of  the  load,  and 
elements  2  and  3  carry  the  other  part.  The  relative  proportions  into  which  Nxy  divides 
are  obtained  by  equating  the  shear  displacement  of  element  1  to  the  sum  of  the  shear 
displacements  of  elements  2  and  3  as  shown  in  figure  3.  The  resulting  internal  shear- 
stress  resultants  are  given  by  equations  (4)  to  (6),  and  the  corresponding  shear  strains 
are  given  by  equation  (7).  For  a  corrugated  panel,  the  stiffness  A00j  •  is  zero  and 

Nxy^  =  Nxy^  =  Nxy. 


During  the  synthesis,  the  internal  stress  resultants  Nx^  and  Nxy^^,  given  by 

equations  (3)  to  (6),  are  compared  with  the  local -buckling  stress  resultants  discussed  in 
a  subsequent  section. 
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Undeformed  position 
Deformed  position 


Figure  3.-  Equating  shear  displacement  of 
element  1  with  sum  of  shear  displacements 
of  elements  2  and  3  to  obtain  internal  shear - 
stress  distribution. 

Stresses  and  strains  in  each  lamina.-  With  0  defined  as  the  angle  between  the  fiber 
direction  and  the  x-direction  (fig.  1),  the  strains  in  each  lamina  in  the  1-direction,  which 
is  parallel  to  the  fiber  direction,  and  in  the  2 -direction,  which  is  perpendicular  to  the  fiber 
direction,  are  given  by  equations  (8)  to  (10).  The  corresponding  stresses. in  each  lamina 
are  given  by  equations  (11)  to  (13). 

The  stresses  and  strains  given  by  equations  (8)  to  (13)  are  used  in  the  strength 
criteria  discussed  in  the  next  section.  Since  each  element  contains  fibers  oriented  in  both 
the  -0  and  -t-6  directions,  stresses  and  strains  at  both  orientations  should  be  examined 
for  strength  considerations  if  shear  loads  exist. 

Strength  Failure  Criteria 

Several  strength  criteria  have  been  proposed  for  designing  filamentary  composite 
laminates.  Two  criteria,  referred  to  here  as  the  maximum-strain  criterion  and  the 
interactive-stress  criterion,  were  used  in  the  studies  included  in  this  report.  (In  no  case 
were  both  criteria  used  simultaneously.)  Both  approaches  appear,  for  example,  in  refer¬ 
ence  8.  The  chosen  failure  criterion  is  applied  to  each  lamina  in  each  element. 

The  strength  failure  criteria  equations  are  summarized  in  table  II  and  are  discussed 
in  the  text. 

Maximum  strain.  -  In  the  maximum -strain  criterion,  upper  and  lower  limits  are 
placed  on  the  permissible  values  of  e  j,  eg,  and  yjg*  upper  limit  is  used  with 
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positive  strains,  and  the  lower  limit  with  negative  strains.  In  terms  of  inequality  expres¬ 
sions,  the  maximum-strain  criterion  is  given  by  equations  (14)  to  (16)  in  which  the  super¬ 
script  a  means  allowable,  and  signs  and  values  of  allowable  strains  must  be  selected  to 
limit  both  positive  and  negative  strains. 

Stress  interaction.  -  While  each  separate  strain  component  is  limited  in  the 
maximum- strain  criterion,  two  of  the  stress  components  are  linked  in  the  stress- 
interaction  criterion  given  here.  Fiber  failure  is  assumed  to  be  prevented  by  equa¬ 
tion  (17),  and  matrix  failure  is  assumed  to  be  prevented  by  equation  (18).  Again,  the 
superscript  a  means  allowable,  and  signs  and  values  of  allowable  stresses  must  be 
selected  to  limit  both  positive  and  negative  stresses. 

Buckling  Analysis 

The  approach  taken  here  is  a  traditional  one  in  which  buckling  is  assumed  to  be 
prevented  if  several  independent  forms  of  buckling,  including  overall-  and  local -buckling 
modes,  are  postulated  and  prevented.  Overall -buckling  equations  are  based  on  smeared 
stiffnesses,  and  local -buckling  equations  are  based  on  simple -support  boundary  conditions 
along  the  lines  of  attachment  of  adjacent  elements.  (This  same  approach  was  used,  for 
example,  in  refs,  9  and  10  and  more  recently  in  ref.  3.)  This  approach  ignores  compati¬ 
bility  of  buckling  mode  shape,  a  factor  that  tends  to  make  the  approach  conservative.  On 
the  other  hand,  the  simplified  analyses  do  not  consider  all  the  potential  buckling  modes,  a 
factor  that  could  make  the  approach  unconservative.  Unpublished  studies  with  a  more 
accurate  buckling  analysis  indicate  that  in  most  cases  the  approach  is  sufficiently  accurate 
to  predict  structural  efficiencies  and  design  trends.  In  some  cases,  such  as  those  in 
which  roll  modes  are  active,  the  buckling  analysis  approach  used  here  is  not  adequate. 

Both  the  overall-  and  local-buckling  analyses  assume  that  the  plate  structure  is 
orthotropic  with  principal  axes  alined  with  the  edges  of  the  plate.  Buckling  caused  by 
longitudinal  compression  and  buckling  caused  by  shear  are  treated  separately.  Combined 
loads  are  accounted  for  with  a  buckling- interaction  formula. 

The  buckling  equations,  which  are  based  on  work  presented  in  references  8  and  11 
to  15,  are  summarized  in  table  III  and  are  discussed  in  the  text.  Certain  detailed  expres¬ 
sions,  including  equations  for  the  smeared  orthotropic  stiffnesses,  are  presented  in  the 
appendix. 

Overall  buckling.  -  The  order  of  presentation  is  longitudinal  compression,  shear,  and 
combined  longitudinal  compression  and  shear. 

Longitudinal  compression!  The  stiffened  panel  is  treated  as  a  wide  column  with  the 
loaded  ends  simply  supported  arid  with  the  lateral  edges  free  as  shown  in  figure  4.  A 
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buckling  equation  that  takes  into  account  transverse  shear  effects  is  presented  as  equa¬ 
tion  (19)  in  table  III.  The  transverse  shearing-stiffness  term  is  taken  as  2  •  b2  •  ^062' 


Figure  4.-  Wide-column  or  Euler  buckling  is  assumed 
for  overall  buckling  in  axial  compression. 


Shear:  The  panel  is  assumed  to  be  long  in  the  y-direction  and  simply  supported 
along  the  edges  x  =  0  and  x  =  L  as  shown  in  figure  5.  The  buckling  equations  are 
given  as  equations  (20)  to  (22)  in  table  III. 

Equations  (21)  and  (22)  are  based  on  the  assumption  that  the  smeared  bending  stiff¬ 
ness  in  the  x-direction  is  much  larger  than  the  smeared  bending  stiffness  in  the 
y-direction  D2,  and  on  the  assumption  that  the  dimension  a  (fig.  5(a))  is  larger  than 
the  dimension  L.  These  assumptions  permit  a  parameter  ^  defined  by 
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1/4 


to  be  taken  as  zero  in  the  derivation  of  equations  (21)  and  (22).  The  equations  correspond, 
therefore,  to  the  infinitely  long  panel  shown  in  figure  5(b). 


5S  =  Simple  support 


(a)  Finite  shear  model.  (b)  Infinitely  long  shear  model. 

Figure  5.  -  Finite  and  infinitely  long  shear  models  for  overall  shear  buckling. 


Combined  longitudinal  compression  and  shear:  When  both  longitudinal  compression 
and  shear  loadings  are  present,  buckling  is  assumed  to  occur  when  the  buckling-interaction 
formula  given  by  equation  (23)  is  satisfied.  Based  on  studies  of  isotropic  plates  in  refer¬ 
ence  14  and  orthotropic  plates  in  reference  15,  this  buckling-interaction  formula  appears 
to  be  adequate  for  the  panels  and  loadings  considered  herein. 

Local  buckling.  -  The  second  level  subscript  (indicating  the  kth  element)  is  dropped 
in  table  III  and  elsewhere  in  this  report  except  where  it  is  necessary  for  clarity. 

Longitudinal  compression:  The  structural  model  is  shown  in  figure  6.  The  length 
L  of  each  element  is  assumed  to  be  much  greater  than  its  width  b.  With  these  simpli¬ 
fying  assumptions  the  buckling  equation,  which  is  applied  to  each  of  the  four  elements  of 
the  cross  section,  is  given  by  equation  (24). 

Shear:  Local  shear  buckling  is  treated  in  the  same  way  that  overall  shear  buckling 
was  treated  in  an  earlier  section.  The  only  difference  is  the  stiffnesses  that  are  used. 
Again,  the  structural  model  is  shown  in  figure  6.  The  equations  for  local  shear  buckling 
are  given  by  equations  (25)  to  (27)  in  table  III. 

Equations  (26)  and  (27)  are  based  on  the  assumption  that  L  »  b  (fig.  6),  This 
assumption  permits  a  parameter  B  defined  by 
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1/4 


L  D 
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to  be  taken  as  zero  in  the  derivation  of  equations  (26)  and  (27), 


Combined  longitudinal  compression  and  shear:  The  buckling- interaction  formula 
given  by  equation  (23)  is  used  for  local  as  well  as  overall  buckling.  The  adequacy  of  that 
formula  is  examined  and  verified  for  orthotropic  plates  in  reference  15. 

Sizing  Technique 

The  design  problem  is  formulated  as  a  nonlinear  mathematical  programing  problem. 
The  design  variables  are  the  element  widths  and  thicknesses  that  define  the  cross  section. 
The  constraints  are  the  strength,  buckling,  and  stiffness  design  requirements.  Since  the 
panel  length  L  is  fixed,  the  objective  function  (the  quantity  that  is  minimized)  is  the 
panel  mass  per  unit  width.  Equations  for  the  strength  and  buckling  constraints  have  been 
presented  in  the  preceding  sections.  Equations  for  the  stiffnesses  are  given  in  the 
appendix. 
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The  specific  design  constraints  considered  herein  are  overall  buckling,  local 
buckling  of  each  element,  strength  of  each  lamina  in  each  element,  extensional  stiffness, 
and  shear  stiffness.  Upper  and  lower  bounds  are  also  placed  on  various  dimensions. 

The  computer  program  used  to  perform  the  optimization  in  the  studies  presented 
herein  is  a  general-purpose  parameter  optimizer  denoted  AESOP  (refs.  16  to  18)  which 
accounts  for  constraints  using  an  exterior  penalty  function  approach.  The  AESOP  code 
contains  several  optimization  algorithms  that  can  be  selected  in  various  combinations  by 
the  user.  The  approach  used  here  was  a  combination  called  "creeper”  and  "pattern." 
These  algorithms  are  described  in  references  16  to  18. 


DESIGN  STUDIES  CARRIED  OUT  WITH  ANALYSIS-DESIGN  PROCEDURE 

Studies  were  carried  out  to  examine  detailed  design  trends,  sensitivity  of  structural 
efficiency  to  changes  in  certain  dimensions,  effect  of  extensional  and  shear-stiffness 
design  requirements,  effect  of  local  buckling  boundary  conditions,  and  the  effect  of  com¬ 
bined  compression  and  shear  loadings.  The  material  properties  that  are  used  for  all 
studies  are  given  in  tables  IV  and  V.  The  upper  and  lower  limits  placed  on  the  cross- 
section  dimensions  are  given  in  table  VI. 


Structural  Efficiency  Diagrams  and  Scaling 


The  structural  design  data  developed  in  this  report  are  presented  in  the  form  of 


structural  efficiency  diagrams  in  which  the  mass  index 


W/A 


N 


xy 


of  minimum-mass  panels 


also  called  the  structural  index. 


Nx 

is  given  as  a  function  of  the  loading  index  or 
Presenting  design  data  in  this  manner  minimizes  size  effects  and  broadens  the  applica¬ 
bility  of  the  data.  If  there  were  no  upper  or  lower  limits  on  the  cross-sectional  dimen¬ 
sions  and  if  the  thickness  could  vary  continuously,  then  size  effects  would  be  eliminated. 
Size  effects  are  eliminated  because  of  dimensional  similarity  or  scaling  principles  that 
apply  to  composite  panels  in  the  same  way  that  they  have  been  applied  to  conventional 
metal  panels.  (See,  for  example,  ref.  19.) 


The  scaling  principle  for  panels  with  loadings  considered  herein  can  be  summarized 
as  follows.  If  two  panels  have  the  same  proportions  -  that  is,  if  all  the  dimensions 
(including  thicknesses  and  length)  of  one  panel  are  multiplied  by  a  single  scale  factor  to 
obtain  the  d^ensions  of  the  other  panel  -  and  if  both  panels  are  subjected  to  loads  which 
give  the  same  value  of  the  loading  index,  then 

(1)  the  stresses  in  the  corresponding  laminas  are  the  same  for  both  panels; 

(2)  the  ratios  of  the  element  stress  resultants  to  the  element  buckling-stress 

resultants  are  the  same  for  both  panels; 


14 


(3)  the  ratio  of  the  overall  stress  resultant  to  the  overall  buckling-stress  resultant 

is  the  same  for  both  panels;  and 

(4)  the  mass  index  is  the  same  for  both  panels. 

These  four  statements  confirm  that  for  a  given  structural  concept  (including  the  material 
system),  minimum-mass  designs  fall  on  a  single  structural  efficiency  curve  and  that 
dimensions  of  minimum-mass  panels  scale  linearly  with  the  load  (not  the  load  index). 

The  results  are  presented  in  terms  of  dimensional  quantities.  In  some  cases,  for 
example,  t2  is  used  as  a  parameter.  In  all  cases  the  data  are  for  panels  with 
L  =  0.76  m  (30  in.).  It  is  understood  that  scaling  techniques  can  be  applied  to  the  results 
to  make  them  applicable  to  other  panel  sizes  and  loads.  For  example,  curves  for  which 

^2 

t2  is  constant  can  be  interpreted  as  curves  for  which  is  constant. 

Hat-Stiffened  Panel,  Longitudinal  Compressive  Loading 
Effect  of  discrete  values  of  thickness  t2(±45°).-  Minimum  values  of  the  mass  index 

'E/A  are  given  in  figure  7  as  a  function  of  compressive  loading  index  ^  for  various 
values  of  t2(±45°).  It  is  assumed  that  the  thickness  of  0°  material  can  vary  continuously. 
The  values  chosen  for  t2  represent  four-ply  and  eight-ply  laminates  made  up  of  plies 
of  either  0.051  mm  (0.0020  in.)  or  0.140  mm  (0.0055  in.)  thickness.  Laminates  that  are 
multiples  of  four  plys  can  be  made  balanced  and  symmetric;  balanced  and  symmetric 
laminates  are  required  in  all  the  studies  presented  in  this  report.  The  curves  give  an 
indication  of  the  range  over  which  a  particular  value  of  t2  is  lightest  and  the  penalty 
involved  in  selecting  another  thickness. 

The  dashed  curve  indicates  designs  for  which  t2(±450)  is  allowed  to  vary  contin¬ 
uously.  That  curve  is  referred  to  as  the  t2free  curve.  It  is  tangent  to  the  various 
curves  for  which  t2  is  constant.  The  designs  associated  with  the  t2free  curve  are 
referred  to  as  the  t2free  designs. 

The  requirement  that  elements  1  and  3  have  a  width  not  less  than  2.03  cm  (0.80  in.) 

N 

begins  to  have  an  effect  for  loadings  less  than  ^  =  276  kPa  (40  Ibf/in^)  where  the 
t2free  curve  begins  to  deviate  from  a  straight  line. 

All  the  examples  presented  herein  use  the  interactive-stress  criterion  given  by 

equations  (17)  and  (18).  For  loadings  greater  than  =  3.45  MPa  (500  Ibf/in^),  strength 
considerations  begin  to  affect  the  t2free  designs  and  the  designs  for  which  t2  =  1.12  mm 
(0.044  in.).  For  this  hat-stiffened  configuration  and  for  the  material  properties  (table  IV) 
used  for  these  studies,  the  maximum  strain  requirement  (eqs.  (14)  to  (16))  produces 
slightly  heavier  designs  than  the  interactive-stress  requirement.  For  example,  at 

—  =  6.89  MPa  (1000  Ibf/in^)  where  the  difference  is  largest,  the  maximum  strain  designs 
L 
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are  about  4  percent  heavier  than  the  interactive-stress  designs.  This  is  the  only  study 
reported  in  this  paper  in  which  the  maximum -strain  criterion  is  used.  Strength  require¬ 
ments  do  not  affect  the  designs  with  the  other  values  of  t2  shown  in  figure  7. 


W/A 

L 


Ibm 


Figure  7.-  Structural  efficiency  of  graphite -epoxy,  hat-stiffened  compres¬ 
sion  panels  for  various  prescribed  values  of  thickness  t2  and  for  t2 
varying  continuously;  L  =  0.76  m  (30  in.);  thickness  t2  is  shown  in 
mm  (in.). 


Consider  the  loading  =  689  kPa  (100  Ibf/in^).  The  lightest  design  has  t2 
equal  to  about  0.41  mm  (0.016  in.).  Manufacturing  requirements  may,  however,  make  it 
more  desirable  to  fabricate  a  design  with  t2  equal  to  0.56  mm  (0.022  in.).  The  cross 
sections  for  the  minimum-mass  designs  having  t2  equal  to  0.41  mm  and  0.56  mm  are 
presented  in  figure  8.  The  design  for  which  t2  =  0.56  mm  is  larger  and  deeper  than 
the  design  for  which  t2  =  0.41  mm.  (Both  designs  are  penalized  slightly  by  the 
b3  =  2.03  cm  (0.80  in.)  requirement.) 
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Period  for  t2  =  0.56  mm  (0.022  in.) - ^ 


— -  ■  t2  =  0.41  mm  (0.016  in.) 

t2  =  0.56  mm  (0.022  in.) 

Figure  8.-  Cross  sections  for  two  hat-stiffened  panels 
N 

designed  for  =  689  kPa  (100  Ibf/in^);  L  =  0.76  m 
(30  in.). 


The  significant  point  here  is  that  even  though  t2  is  increased  from  0.41  mm 
(0.016  in.)  to  0.56  mm  (0.022  in.)  -  a  37-percent  increase  -  the  design  for  which 
t2  =  0.56  mm  is  only  4 -percent  heavier  than  the  design  for  which  t2  =  0.41  mm.  In  this 
case,  as  well  as  in  many  other  cases,  the  manufacturing  requirement  on  the  thickness  t2 
can  be  accommodated  with  very  little  penalty  in  structural  efficiency.  It  should  be  pointed 
out,  however,  that  although  there  is  very  little  difference  in  the  mass  of  the  two  designs, 
the  extensional  and  the  shearing  stiffnesses  of  the  one  design  differ  substantially  from  the 
corresponding  stiffnesses  in  the  other  design.  This  fact  is  discussed  in  more  detail  in  a 
subsequent  section. 


Variation  of  dimensions  as  function  of 


5c 

L 


for 


t2  =  0.56  mm  (0.022  in.).-  Cross- 


section  dimensions  of  the  designs  in  figure  7  for  which  t2  =  0.56  mm  (0.022  in.)  are 

Nx 

presented  as  a  function  of  the  loading  index  —  in  figure  9.  These  dimensions  are  for 
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Figure  9.-  Variation  of  cross-sectional  dimensions  as  function  of 
loading  for  minimum-mass,  hat -stiffened,  graphite-epoxy 
compression  panels;  tn  =  0.56  mm  (0.022  in.);  L  =  0.76  m 


a  panel  with  a  length  L  of  0.762  m  (30  in.)  and  for  the  material  properties  of  table  IV. 
The  loading  ranges  over  which  various  buckling  constraints  are  active  are  given  at  the  top 
of  figure  9.  Element  3  is  not  buckling  critical  over  the  loading  range  shown,  in  figure  9. 

The  width  b3  of  element  3  is  shown  in  figure  9  as  being  equal  to  its  minimum 

N 

value  (2.03  cm  (0.80  in.))  until  is  equal  to  about  4830  kPa  (700  Ibf/in^).  Actually, 

^  Nx 

b3  has  a  weak  tendency  to  increase  to  a  larger  value  over  the  range  =  2410  to 
3450  kPa  (350  to  500  Ibf/in^).  This  increase  in  b3  is  accompanied  by  a  decrease  in 
t3(0O)  so  that  at  a  given  loading  the  total  longitudinal  extensional  stiffness  (ET3  •  b3)  of 
element  3  remains  almost  constant  regardless  of  the  relative  values  of  b3  and  13(0°). 
The  simple  variation  of  b3  and  13(0°)  is  shown  in  figure  9  for  clarity.  The  weak 
tendency  of  b3  to  exceed  its  minimum  value  in  the  above  loading  range  is  discussed  in 
the  next  section. 

Cross-section  shapes  of  the  panels  having  13  =  0.56  mm  (0.022  in.)  and  designed  to 
N 

support  loadings  of  =  68.9,  689,  and  6890  kPa  (10,  100,  and  1000  Ibf/in^)  are  pre- 

sented  in  figure  10.  The  designs  for  =  68.9  and  6890  kPa  are  shown  merely  to 
indicate  the  design  trends  for  a  fixed  value  of  t2.  Other  values  of  t2  provide  substan¬ 
tially  lighter  designs  at  these  two  loadings. 

,  ^  0  1  2  3 

Inches  1 _ [ _ I _ } 

0  2  4  6  8 

Centimeters  I _ I — I — I — i — I — I — 1 — 1 


"^7  ^ 

(a)  ^  =  68.9  kPa  (10  Ibf/in^). 

JL 


(b)  =  689  kPa  (100  lbf/in2). 

L 


(c)  ^  =  6890  kPa  (1000  Ibf  (in^). 

Figure  10.-  Cross-sectional  shapes  of  hat- stiffened,  graphite -epoxy  compression 
panels  having  t2  =  0.56  mm  (0.022  in.)  and  L  =  0.76  m  (30  in.)  for  three 
values  of  axial  compressive  load. 
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Effect  of  changes  in  width  of  element  3.  -  In  the  previous  section  it  was  pointed  out 
that  although  the  width  b,  of  element  3  is  shown  in  figure  9  as  being  equal  to  the  mini- 

mum  value  (2.03  cm  (0.80  in.))  until  ^  is  equal  to  about  4830  kPa  (700  Ibf/in^),  there 

Li 

is  actually  a  weak  tendency  for  b3  to  be  larger  than  the  minimum  value  in  the  region 
N 

=  2410  to  3450  kPa  (350  to  500  Ibf/in^).  That  effect  is  shown  in  figure  11  for  a  loading 

Li 

of  ^  =  2760  kPa  (400  Ibf/in^).  In  figure  11  stiffened  panels  are  completely  redesigned 
for  various  prescribed  values  of  b3;  that  is,  each  point  shown  by  a  circular  symbol 
represents  a  minimum-mass  design  for  that  value  of  b3.  For  these  calculations, 
t2  =  0.56  mm  (0.022  in.).  The  minimum  value  of  the  mass  index  —  -  occurs  at  b3 
equal  to  about  2.54  cm  (1.0  in.).  For  this  loading,  the  penalty  for  choosing  bg  =  2.03  cm 

Nv 

is  about  0.25  percent.  It  is  apparent  that  for  this  value  of  — ,  the  width  b3  can  vary 
over  a  large  range  with  very  little  change  in  structural  efficiency.  The  width  bg  is  a 
weak  design  variable  for  other  loadings  also.  This  fact  can  be  exploited  to  accommodate 
manufacturing  or  other  design  considerations. 


Percent  change  in 


W/A 


100- 


W/A  _  /W/A\ 

L  V  L  /min 

(m\ 

\  L  /  min 


I _ I _ 1 _ I _ I _ I 

2.0  2.5  3.0  3.5  4.0  4.5  cm 

I _ I _ ^ ^ _ \ _ I 

.8  1.0  1.2  1.4  1.6  1.8  in. 

^3 

Figure  11.-  Percent  change  in  structural  efficiency  as  function  of 
dimension  bg  for  minimum-mass,  graphite -epoxy,  hat-stiffened 
compression  panels;  t2  =  0.56  mm  (0.022  in.);  L  =  0.76  m 

(30  in.);  ^  =  2760  kPa  (400  Ibf/in^). 


20 


Effect  of  changes  in  width  of  element  4.  -  The  structural  efficiency  is  more  sensitive 
to  changes  in  than  to  changes  in  bg.  An  indication  of  that  sensitivity  is  shown  in 
figure  12  in  which  the  percent  change  in  mass  index  is  presented  as  a  function  of  width  d 

(twice  of  the  flat  portion  between  stiffeners  for  two  values  of  689  and  4830  kPa 

(100  and  700  Ibf/in^).  The  circular  and  square  symbols  indicate  minimum-mass  designs 
for  the  689  and  4830  kPa  cases,  respectively.  For  these  results,  the  thickness  t2  is 
allowed  to  vary  freely. 


d,  in. 


d,  cm 


Figure  12.  -  Effect  on  structural  efficiency  of  increasing  stiffener 
spacing  d  of  hat -stiffened,  graphite -epoxy  compression  panels; 
t2  varies  continuously;  L  =  0,76  m  (30  in,). 

First,  consider  the  case  for  which  =  689  kPa  (100  Ibf/in^).  With  no  constraint 
placed  on  d,  the  lightest  design  occurs  for  d  equal  to  about  4,8  cm  (1,9  in,).  As  d 
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W /A 

is  increased,  the  mass  index  — increases.  When  d  is  equal  to  about  7.9  cm 
(3.1  in.),  the  lightest  design  changes  from  one  with  0^  material  in  elements  3  and  4  to  one 
that  has  0^  material  in  element  3  only  ^t^(O^)  =  0^.  The  two  concepts,  one  with  and  the 
other  without  0^  material  in  element  4,  do  not  evolve  from  one  another.  They  are  com¬ 
pletely  independent.  The  concept  with  no  0^  material  in  element  4  is  similar  to  those  for 
low  loadings  shown  in  figure  9  in  that  bj  is  equal  to  b4.  As  d  is  increased  from 
about  4.8  cm  to  12.7  cm  (5.0  in.),  the  mass  index  increases  by  about  42  percent. 

For  the  case  in  which  =  4830  kPa  (700  Ibf/in^),  the  efficiency  of  the  concept  for 
which  14(0^)  =0  is  restricted  because  of  the  upper  bound  placed  on  tg(O^)  and  because 
of  strength  considerations.  For  these  reasons,  the  concept  with  0^  material  in  both 
elements  3  and  4  is  lighter  than  the  concept  with  0^  material  in  only  element  3  for  the 
entire  range  of  d  considered.  The  lightest  design  occurs  for  d  equal  to  about  8. 1  cm 
(3.2  in.).  As  d  is  increased  to  12.7  cm  (5.0  in.),  the  mass  index  increases  by  about 
18  percent. 

Another  way  of  studying  the  effect  of  increasing  the  distance  d  between  stiffeners 
is  with  a  structural  efficiency  diagram  such  as  that  shown  in  figure  13.  The  objective 
here  is  to  show  the  effects  of  a  constraint  defined  by  d  ^  12.7  cm  (5.0  in.).  Since  there 
is  no  tendency  for  d  to  be  greater  than  12.7  cm  for  the  cases  considered,  the  constraint 
becomes  simply  d  =  12.7  cm. 

The  bottom  solid  curve  in  figure  13  gives  the  structural  efficiency  of  hat- stiffened 
designs  with  no  constraint  on  d.  It  is  the  t2free  curve  of  figure  8.  The  two  dashed 
curves  that  are  tangent  to  the  t2free  curve  are  portions  of  the  corresponding  curves 
of  figure  8.  Like  the  solid  curve,  these  two  dashed  curves  have  no  constraint  on  d  and 
are  used  for  reference.  For  the  other  curves,  d  is  required  to  be  12.7  cm  (5.0  in.). 

The  other  solid  curve  in  figure  13  forms  the  lower  bound  for  the  designs  having 
d  =  12.7  cm  (5,0  in.).  For  this  curve,  t2  is  allowed  to  vary  freely.  For  loadings  less 

than  about  =  3450  kPa  (500  Ibf/in^),  the  designs  for  which  t4(0^)  =  0  are  lighter  than 
the  designs  for  which  14(0^)  0.  However,  because  of  the  upper  bound  on  t3(0^)  and 

because  of  strength  considerations,  the  designs  for  which  14(0^)  ^  0  are  lighter  for 

N 

loadings  greater  than  about  =  3450  kPa  (500  Ibf/in^). 

L 

The  penalty  in  structural  efficiency  caused  by  increasing  the  distance  d  between 

N 

stiffeners  to  12.7  cm  (5.0  in.)  depends  upon  the  loading  level.  At  =  68.9  kPa 

L 

,  ,  9 

(10  Ibf/ln'^),  the  increase  in  the  mass  index  is  about  77  percent;  at  —  =  6890  kPa 
(1000  Ibf/in^),  the  increase  is  about  14  percent. 
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Figure  13.-  Effect  of  d  =  12.7  cm  (5  in.)  stiffener  spacing  requirement  on 
structural  efficiency  of  hat-stiffened,  graphite -epoxy  compression  panels. 
Thicknesses  are  shown  in  mm  (in.);  L  =  0.76  m  (30  in.). 

Effect  of  local -buckling  boundary  conditions.  -  As  explained  earlier,  the  analysis  for 
local  buckling  is  based  on  the  assumption  that  each  element  is  simply  supported  along  its 
line  of  attachment  to  the  adjacent  elements.  This  assumption  is  believed  to  be  conserva¬ 
tive  in  most  cases.  In  particular,  since  element  3  never  experiences  local  buckling  in  the 
designs  generated  in  this  report,  it  is  reasonable  to  assume  that  element  3  provides  rota¬ 
tional  support  to  element  2  and  that  a  more  sophisticated  analysis -design  procedure  could 
exploit  that  extra  support 

To  assess  the  effect  of  local -buckling  boundary  conditions  and  to  provide  information 
on  the  potential  benefits  available  with  a  more  sophisticated  analysis,  additional  designs 
were  generated  under  the  assumption  that  the  local -buckling  boundary  conditions  for 
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elements  1  and  2  were  clamped  rather  than  simply  supported.  The  results  are  presented 
in  figure  14.  The  solid  curves  represent  designs  obtained  with  the  clamped  boundary 
condition  assumption,  and  the  dashed  curves  are  the  original  designs  with  the  simple- 
support  boundary  conditions. 


Figure  14.  -  Effect  of  local-buckling  boundary  conditions  on  structural  efficiency 
of  hat -stiffened,  graphite -epoxy  compression  panels;  L  =  0.76  m  (30  in.). 
Thicknesses  are  shown  in  mm  (in.). 

The  benefits  derivable  from  the  clamped  boundary  condition  depend  upon  the  curves 
being  compared  and  the  loadings  at  which  the  curves  are  compared.  In  the  loading  region 
where  the  designs  with  t2  =  0.56  mm  (0.022  in.)  are  the  most  efficient  designs,  clamping 
provides  a  10-percent  reduction  in  the  mass  index.  These  results  suggest  that  a  more 
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sophisticated  buckling  analysis,  such  as  those  described  in  references  20  and  21,  would 
provide  designs  that  are  about  10  percent  more  efficient  than  those  presented  in  this 
report. 

Stiffness  considerations.  -  It  was  shown  in  reference  5  that  the  longitudinal  exten- 
sional  stiffness  of  minimum-mass  graphite -epoxy  compression  panels  can  be  less  than 
that  of  minimum-mass  aluminum  panels  designed  to  carry  the  same  load.  This  fact  is 
demonstrated  in  figure  15  in  which  the  longitudinal  extensional  stiffness  ET  is  presented 

Nx 

as  a  function  of  the  loading  index  -j—  for  a  portion  of  the  designs  shown  in  figure  7. 


Figure  15.-  Longitudinal  extensional  stiffness  as  function  of  loading  for 
minimum-mass,  graphite -epoxy  and  aluminum  hat-stiffened  compression 
panels;  L  =  0.76  m  (30  in.).  Thicknesses  are  shown  in  mm  (in,). 


The  solid  lines  represent  graphite -epoxy  panels  with  four  values  of  t2,  and  the  dashed 
line  represents  aluminum  panels.  The  discontinuities  in  the  extensional  stiffness  of  the 
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graphite -epoxy  designs  occur  where  the  most  efficient  designs  shift  from  one  value  of  t2 
to  another.  This  relationship  can  be  seen  by  comparing  figure  15  with  figure  7.  The  main 
reason  for  the  "sawtooth"  effect  in  figure  15  is  the  variation  of  0°  material  in  the  panels. 

If  stiffness  is  an  important  design  consideration,  that  fact  can  be  taken  into  account  by 
including  stiffness  constraints  in  the  synthesis  procedure. 

Aluminum  wing  panels  on  commercial  jet  transports  may  have  a  much  higher  longi¬ 
tudinal  extensional  stiffness  than  that  shown  in  figure  15,  (See  ref,  6,)  Commercial  air¬ 
craft  may  also  have  large  shear -stiffness  requirements  as  well  as  stiffener  spacing 
requirements.  To  study  the  effect  of  these  additional  design  requirements,  minimum- 
mass  panels  were  designed  for  various  combinations  of  extensional  and  shear  stiffness. 
(Expressions  for  both  stiffnesses  are  given  in  the  appendix,)  All  panels  were  designed  to 

carry  a  load  of  ^  =  689  kPa  (100  Ibf/in^),  and  all  panels  have  a  distance  d  between 
stiffeners  of  12.7  cm  (5.0  in.)  and  a  length  L  of  0.76  m  (30  in.).  The  thickness  tg  is 
allowed  to  vary  freely. 

W/A 

The  results  are  presented  in  figure  16  in  which  the  mass  index  of  the 

minimum-mass  panels  is  shown  as  a  function  of  the  design  requirements  on  shear  stiff¬ 
ness  and  longitudinal  extensional  stiffness.  The  dashed  lines  are  lines  of  constant  value 
of  the  mass  index.  The  circular  symbol  represents  the  design  which  has  no  stiffness 
requirement.  The  solid  line  above  the  symbol  represents  the  locus  of  minimum-mass 
designs  that  meet  an  extensional-stiffness  requirement,  and  the  solid  line  below  and  to  the 
right  of  the  solid  symbol  represents  the  locus  of  minimum-mass  designs  that  meet  a 
shear -stiffness  requirement.  Points  within  the  region  containing  the  dashed  lines  repre¬ 
sent  designs  with  combined  extensional-  and  shear-stiffness  requirements.  The  discon¬ 
tinuity  in  the  line  defining  the  extensional-stiffness  designs  is  caused  by  a  change  in 
design  concept  similar  to  that  shown  in  figure  12. 

Two  observations  follow  from  the  results  presented  in  figure  16,  First,  starting 
at  the  symbol,  increasing  either  stiffness  with  a  minimum  increase  in  mass  (that  is, 
moving  along  the  solid  line)  causes  the  other  stiffness  to  decrease.  Except  for  the  lower 
boundary,  this  characteristic  exists  for  all  designs.  This  characteristic  exists  primarily 
because  the  extensional  stiffness  is  increased  by  adding  0®  material  and  the  shear  stiff¬ 
ness  is  increased  with  ±45°  material.  Second,  the  penalty  in  structural  efficiency 
associated  with  increasing  the  shear  stiffness  is  much  greater  than  the  penalty  associated 
with  increasing  the  extensional  stiffness. 
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Figure  16.-  Panel  structural  efficiency  as  function  of  longitudinal-extensional- 
stiffness  and  shear-stiffness  requirements  for  minimum-mass,  hat -stiffened, 
graphite -epoxy  compression  panels  designed  to  support  loading  of 

=  689  kPa  (100  Ibf/in^).  Thickness  t2  varies  continuously.  Mass 
index  — given  in  kg/m^  (Ibm/in^). 


Hat-Stiffened  Panel,  Combined  Shear  and  Compression 

The  effect  of  a  shear  loading  on  the  structural  efficiency  of  hat-stiffened  compres- 

W  /A 

Sion  panels  is  shown  in  figures  17  and  18.  The  mass  index  — y —  is  presented  as  a 

^x  Nxv  ^  Nxv 

function  of  -rp—  for  several  values  of  in  figure  17  and  for  several  values  of 

L  L  Nx 

in  figure  18.  For  each  curve,  t2  is  allowed  to  vary  continuously,  so.  the  cusps  associated 

with  constant  values  of  t2  do  not  appear.  As  is  to  be  expected,  the  larger  the  ratio  of 

^xv  ^x 

to  — ,  the  larger  the  effect  on  the  mass  index. 

Li  L 
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Figure  18.-  Structural  efficiency  of  graphite -epoxy,  hat-stiffened  panels 
subjected  to  combinations  of  longitudinal  compression  and  shear  loadings; 
L  =  0.76  m  (30  in.);  to  varies  continuously. 


Corrugated  Panel,  Longitudinal  Compressive  Loading 

To  study  the  effect  of  discrete  thickness,  corrugated  compression  panels  were 
designed  with  various  thicknesses  of  ±45°  material.  The  results,  which  are  presented  in 
figure  19,  show  the  same  general  trends  as  the  hat- stiffened  panels  (fig.  7).  The  two 
main  differences  are:  first,  the  values  of  the  mass  index  are  lower  for  the  corrugated 
panel  than  the  hat-stiffened  panel;  second,  the  most  efficient  loading  range  for  a  given 
value  of  t2  is  lower  for  a  corrugated  panel  than  for  a  hat-stiffened  panel. 


-p.kPa 

Figure  19.-  Structural  efficiency  of  corrugated,  graphite -epoxy  compression 
panels  for  various  values  of  thickness  t2;  L  =  0.76  m  (30  in.). 
Thicknesses  are  shown  in  mm  (in.). 
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Corrugated  Panel,  Combined  Shear  and  Compression 

The  effect  of  a  shear  loading  on  the  structural  efficiency  of  corrugated  compression 
panels  is  shown  in  figures  20  and  21.  The  curves  are  similar  to  those  presented  in 
figures  17  and  18.  The  main  difference  is  that  the  mass  index  is  lower  for  the  corrugated 
panels  than  for  the  hat- stiffened  panels. 
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Figure  20.-  Effect  of  shear  loading  on  structural  efficiency  of  corrugated, 
graphite -epoxy  compression  panels;  L  =  0. 76  m  (30  in.);  t2  varies 
continuously. 
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Figure  21.-  Structural  efficiency  of  graphite -epoxy  corrugated  panels 
subjected  to  combinations  of  longitudinal  compression  and  shear 
loadings;  L  =  0.76  m  (30  in.);  t2  varies  continuously. 

Hat-Stiffened  and  Corrugated  Panels,  Shear  Loading 

The  structural  efficiency  of  hat-stiffened  and  corrugated  panels  loaded  only  by  a 
shear  loading  is  shown  in  figure  22.  The  solid  curves  are  for  graphite -epoxy  panels. 

Over  the  range  of  loading  considered,  the  graphite -epoxy  corrugated  panels  are  5  to 
20  percent  lighter  than  the  graphite -epoxy  hat- stiffened  panels. 

The  dashed  curve  is  for  aluminum  corrugated  panels  having  uniform  thickness  and 

designed  using  the  same  procedure  as  that  used  for  the  composite  panels.  The  properties 

assumed  for  the  aluminum  material  are  given  in  table  V.  In  the  lower  loading  range 
N 

<  345  kPa  (50  Ibf/in^),  where  the  panels  are  buckling  critical  but  not  strength  critical, 
the  aluminum  corrugated  panels  are  about  45  percent  heavier  than  the  graphite -epoxy 
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corrugated  panels.  For  loadings  greater  than  =  345  kPa,  material  strength  con- 

Li 

siderations  become  dominant  for  the  aluminum  panels,  and  these  panels  become  consider¬ 
ably  heavier  than  the  graphite -epoxy  panels. 


kPa 

Figure  22.-  Effect  of  shear  loading  on  structural  efficiency  of  graphite -epoxy 
and  aluminum  corrugated  panels  and  graphite -epoxy  hat- stiffened  panels; 
L  =  0. 76  m  (30  in.);  t2  varies  continuously. 


The  symbol  in  figure  22  represents  the  theoretical  structural  efficiency  of  an 
aluminum  corrugated  shear  panel  given  in  reference  22.  The  structural  efficiency  of  that 

panel  is  shown  by  the  intersection  of  the  dashed  curves  (overall-  and  local -buckling 

s  /  s 

curves)  at  -^  =  0.05  ksi  in  figure  9  of  reference  22.  (The  quantity  is  the  same  as 

XT  \  u/2  \ 


XV  \  \ 

-^.1  The  structural  efficiency  is  influenced  by  a  12 -percent  increase  in  mass  caused  by 


doublers  at  the  edges  x  =  0  and  x  =  L  and  by  edge  constraints  that  provide  more 
support  than  simple -support  boundary  conditions.  That  panel  is  about  3  percent  heavier 
than  the  aluminum  corrugated  panel  designed  with  the  procedure  presented  in  this  report. 

CONCLUDING  REMARKS 

An  engineering  design  procedure  for  sizing  composite  hat-stiffened  and  corrugated 
panels  subjected  to  longitudinal  compression  and  shear  loadings  has  been  developed  and 
exercised.  The  procedure  is  based  on  nonlinear  mathematical  programing  techniques. 
Expressions  are  given  for  the  design  constraints  used  in  the  procedure.  These  constraints 
are:  overall-buckling  loads,  local -buckling  loads,  lamina  strength,  longitudinal  exten- 
sional  stiffness,  and  shear  stiffness. 

The  buckling  analysis  is  simplified  in  that  only  certain  specific  buckling  modes  are 
considered.  Because  explicit  eigenvalue  solutions  are  available  for  these  buckling  modes, 
eigenvalue  analyses  are  not  required.  Design  studies  covering  a  broad  range  of  design 
parameters  can,  therefore,  be  carried  out  with  a  relatively  small  amount  of  computer 
time.  Unpublished  studies  with  a  more  accurate  buckling  analysis  indicate  that  in  most 
cases  the  approach  is  sufficiently  accurate  to  predict  structural  efficiencies  and  design 
trends.  In  some  cases,  such  as  those  in  which  roll  modes  are  active,  the  buckling  analysis 
approach  used  here  is  not  adequate. 

In  order  to  obtain  a  clearer  understanding  of  both  the  analysis-design  procedure  and 
the  composite  panel  design  results  that  could  be  obtained  by  use  of  that  procedure,  several 
design  studies  were  carried  out.  Buckling  and  strength  design  requirements  are  included 
in  all  the  studies.  In  one  study  the  effect  of  stiffness  requirements  is  explored.  Struc¬ 
tural  efficiency  diagrams  are  used  to  present  the  mass  for  panels  designed  to  support 
longitudinal  compressive  loads,  shear  loads,  and  the  increase  in  mass  required  for 
combined  shear  and  compressive  loadings.  T-he  studies  indicate  that  panel  mass  is  not 
extremely  sensitive  to  panel  dimensions;  that  is,  panels  with  dimensions  moderately 
different  from  those  of  the  lightest  panels  may  be  only  slightly  heavier  than  the  lightest 
panels.  This  fact  can  be  exploited  to  accommodate  manufacturing  and  other  design  con¬ 
siderations.  It  was  also  determined  that  scaling  techniques  in  common  use  for  sizing 
conventional  metal  panels  apply  equally  well  for  composite  panels  where  the  number  of 
design  variables  is  large. 

Other  observations,  based  on  limited  studies,  are: 

1.  There  is  a  greater  mass  penalty  associated  with  increasing  the  shear  stiffness  of 
hat-stiffened  panels  than  with  increasing  the  extensional  stiffness. 
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2.  The  mass  of  hat- stiffened  panels  is  affected  more  by  changes  in  the  distance 
between  stiffeners  than  in  the  width  of  the  hat  cap. 

3.  As  has  been  found  in  studies  for  compression  panels,  composite  shear  panels 
offer  substantial  mass  savings  over  metal  shear  panels. 

The  design  procedure  presented  in  this  report  should  be  used  in  conjunction  with 
engineering- judgment  and  other  analyses  to  insure  the  adequacy  of  the  final  design.  Other 
design  requirements  arising  from  manufacturing  considerations,  for  example,  could  also 
have  an  impact  on  the  design.  To  the  fullest  extent  possible,  all  quantifiable  requirements 
should  be  considered  early  in  the  synthesis. 

Langley  Research  Center 

National  Aeronautics  and  Space  Administration 
Hampton,  VA  23665 
July  23,  1976 
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APPENDIX 


STIFFNESS  EXPRESSIONS  APPEARING  IN 
ANALYSIS-DESIGN  PROCEDURE 


Basic  Laminate  Constitutive  Relations 
The  plate  constitutive  equation  can  be  written  as 


N, 


Ai2 

Ai6 

Bn 

Bi2 

Bie 

A22 

A26 

Bi2 

B22 

B26 

A26 

A66 

B16 

B26 

Bee 

^12 

B16 

Dll 

Di2 

Die 

^22 

^26 

Di2 

D22 

D2e 

^26 

^66 

Die 

D26 

Dee 

where  Ajj ,  Btj ,  and  Djj  are  laminate  stiffnesses  that  are  calculated  from  basic  mate¬ 
rial  properties  by  using  standard  expressions  found,  for  example,  in  reference  8.  The 
quantities  Mx,  My,  and  Mxy  are  plate  moments,  and  kx,  ky,  and  kxy  are  the  asso¬ 
ciated  curvatures.  In  all  cases  considered  in  this  report,  the  ply-stacking  sequence  is 
balanced  and  symmetric  so  that  the  A10,  A20,  and  Bjj  terms  are  zero.  Also,  the 
Dj^0  and  D20  terms  are  assumed  to  be  zero.  For  these  types  of  laminates,  which  are 
denoted  "specially  orthotropic,"  the  set  of  equations  represented  by  equation  (Al)  uncouple 
and  reduce  to 


Nx  All  Ai2  0  €x 


Ny  =  Ai2  A22  0  fiy 


0  0  A00  r. 
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Mx  1^11  ^12  ® 

My  =  D]^2  ^22  ® 


0  0  Dgg  kxy 


In  subsequent  equations,  the  notations  Aij^^  and  Dijj^  denote  the  values  of  Ay  and 


Dy,  respectively,  for  element  k. 


Extensional  Stiffness 


If  the  transverse  load  Ny  is  zero,  the  longitudinal  extensional  stiffness  ETj^  of 
element  k  is  given  by 

If  it  is  assumed  that  all  elements  of  the  panel  cross  section  have  the  same  longitudinal 
strain  ex,  the  longitudinal  extensional  stiffness  ET  is  given  by 

g-p  _  ®^lbi  2ET2b2  +  ET3b3  +  2ET4b4 
bj  +  2b4 

For  the  corrugated  panel,  ETj  is  zero. 

Shear  Stiffness 


The  shear  stiffness  GT  of  the  panel  is  defined  as 


GT  = 


in  which  is  the  average  shear  angle  given  by 


^yxy2^2  ^xys^S  ^l^xy4^4 


bj  +  2b4 


With  equations  (4)  to  (7)  used  to  calculate  the  individual  shear  angles,  the  shear  stiffness 
GT  is 


'2b2  ^  b3  ' 
^602  -^663 


b  j  +  2b^ 


^'^2^64  ^3^664  •'^664 


^1^662 
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Equation  (A8)  can  be  used  for  either  hat-stiffened  or  corrugated  panels.  For  hat-stiffened 
panels  only,  the  following  equation  is  equivalent  to  equation  (A8)  and  provides  more  insight 
into  the  effect  of  the  hat  on  the  stiffness  GT: 


GT  = 


bj  +  2b4 


2b2A66j  b3Ae6j\  Ae64 

^^1^662  ^1^663/ 


Smeared  Stiffnesses  for  Overall  Buckling 

The  overall-buckling  equations  given  in  table  III  are  based  on  the  assumption  that 
the  stiffened  panel  is  an  orthotropic  plate  with  smeared  stiffness  properties.  The  smeared 
stiffnesses  are  defined  in  the  following  differential  equation  for  lateral  deflections  of  an 
orthotropic  plate  loaded  as  shown  in  figure  1: 

. 2/d^,  +  ■■  +  — y—  ^ 


PxMy  9x 


MxMy/  ax^ay 


1 


^bcMy  ay^ 


“^^xy  ax  ay 


(AlO) 


The  stiffnesses  that  appear  in  equation  (AlO)  are  related  to  the  stiffnesses  that  appear  in 
equations  (20)  to  (23)  by 


(All) 


1  -  MxMy 


(A12) 


“3  -  "V  *  jTin 


(A13) 


The  objective  of  this  section  is  to  provide  expressions  for  El,  Dj,  03,  and  D3. 
The  approach  is  the  same  as  that  used  in  reference  23. 

Calculations  presented  in  reference  23  show  that  for  metal  corrugation- stiffened 
panels,  the  product  /%My  can  be  neglected  when  compared  with  1.0.  It  is  assumed  that 
the  product  can  also  be  neglected  for  composite  panels  of  the  type  considered  in  this 
report. 
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Bending  stiffnesses  El  and  Dj.-  The  panel  bending  stiffness  is  simply  the 

El  per  unit  width  of  a  beam  having  the  cross-sectional  shape  of  the  panel.  Since  Pxi^y 
is  small  compared  with  1.0,  Dj  is  equal  to  D^.  It  is  first  necessary  to  locate  the  neu¬ 
tral  axis  of  the  cross  section.  The  distance  from  the  center  of  element  3  to  the  neutral 
axis  of  element  k  is  denoted  zj^.  These  distances  are  shown  in  figure  23  and  are  sum¬ 
marized  below. 


Figure  23.-  Additional  dimensions  for  hat-stiffened  panel. 


Hat-Stiffened  Panel 

Corrugated  Panel 

=  D  +  t^(O^)  +  t2(±45^) 

=  D 

(A14) 

Z2  =  D/2 

Z2  =  D/2 

(A15) 

Zg  =  0 

Zg  =  0 

(A16) 

z^  =  D  +  ^jt4(0°)  +  t2(±45®^ 

Z4  =  D 

(A17) 

The  distance  z  from  the  center  of  element  3  to  the  neutral  axis  of  the  panel  cross  sec¬ 
tion  is  then  given  by 


z  = 


(A18) 
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in  which  ETj^  is  given  by  equation  (A4).  In  equation  (A18),  proper  account  must  be  taken 
of  the  fact  that  elements  2  and  4  appear  twice  in  one  period,  and  the  summations  are  over 
one  period.  The  bending  stiffness  El  of  one  period  of  the  panel  cross  section  is 

The  bending  stiffness  of  each  element  about  its  own  centroid  is  neglected  in  equation  (A19). 
The  error  introduced  by  this  assumption  is  only  about  1  percent. 

The  smeared  bending  stiffnesses  and  Dj  are  the  El  per  unit  width  of  the 

panel 

Dl  -  Dx  -  ^  (A20) 

in  which  p  is  the  period  bj  +  2b^. 

Bending  stiffness  D2.-  The  approach  presented  in  reference  23  is  again  followed. 

The  smeared  bending  stiffness  in  the  y-direction  is  given  by 


where 


and 


(A21) 


(A22) 


4 

sin  a 

b3 

D 

2 

^222 

2 

^223 

sin^a 

^^3 

3 

^222 

2 

^223 

(A23) 


Angle  a  is  the  angle  that  element  2  makes  with  the  horizontal  as  shown  in  figure  23. 

Stiffness  D3.-  For  the  hat-stiffened  panel,  MyD^  can  be  neglected  when  compared 
with  Dxy.  In  this  case  D3  =  D^y.  The  twisting  stiffness  Djjy  is  calculated  from  Bredt 
theory  (see,  for  example,  ref.  24)  and  is  given  by 
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^3  "  ^xy  - 


2a2 


bi  2bo 
p/— +  + 


^66^  ^002  A00g 


in  which  A  is  the  area  enclosed  by  the  hat  and  is  given  by 
A  =  i(bi  +  b3)[D  +  t4(0O)  +  t2(±45°)] 


(A24) 


(A25) 


Experimental  work  presented  in  reference  23  and  analytical  studies  carried  out  to 
support  reference  25  indicate  that  the  value  of  D^y  given  in  equation  (A24)  is  too  large. 
In  the  studies  presented  in  this  report  the  twisting  stiffness  was  obtained  by  applying  a 
0.3  factor  to  the  value  given  by  equation  (A24). 

In  the  case  of  the  corrugated  panel,  /iyE>x  cannot  be  ignored  when  compared  with 
Dxy.  The  stiffness  D3  is  given  by 


^  _  2b2(Di22  +  2000^)  +  2b3(bi22  +  2000^) 


(A26) 
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TABLE  I.-  EQUATIONS  FOR  STRESS  ANALYSIS 


Loading-equation 

identification 

Longitudinal 

compression 

Shear 

Lamina  strains 

referred  to 
lamina  principal 

axes 

y 

Lamina  stresses 

referred  to 
lamina  principal 

Equation 


Equation 

number 


= 


Nx 

ET 


12, 


e,,  =  -e 


yk  ^^A 


22, 


Nxk  "  ^x®\ 


(1) 

(2) 

(3) 


^xy4  -  ^xy 


N_,.  =  N, 


N 


xy 


^2  ^'^3  2b2A66i  ^  bgAeei 

*^1  ^662  ”1^663 


^xy^  “  ^xy  “  Nxyg 


N 


y 


^yk  A66k 


(4) 

(5) 

(6) 

(7) 


€  J  =  cos2  0  ex  +  Sin2  0  ey  +  cos  0  sin  6 
Sg  =  sin^  ^  +  cos^  0  €y  -  cos  0  sin  0  y 


xy 


=  -2  sin  0  cos  0  ex  +  2  sin  0  cos  0  6y 
+  (cos^  0  -  sin3  0)  y^ 


(8) 

(9) 

(10) 


O'!  =  Qii^i  +  Qi2^2 

^2  ^  *^12^1  ^22^2 


(11) 

(12) 


axes 


'’'12  "  12 


(13) 


TABLE  n.-  EQUATIONS  FOR  STRENGTH  FAILURE  CRITERIA 

Failure 

criterion 

Equation 

(*) 

Equation 

number 

Maximum 

strain 

ilsi 

(14) 

2 

(15) 

^<1 

ya 

12 

(16) 

Stress 

interaction 

4  =  1 
^1 

(17) 

/  x2  \2 

/!2\  Jll2  U  1 

La  ]  - 

V^2/  \  12/ 

(18) 

Superscript  a  denotes  allowable  value, 


TABLE  m.-  EQUATIONS  FOR  OVERALL  AND  LOCAL  BUCKLING 


TABLE  IV,-  PROPERTIES  OF  GRAPHITE -EPOXY  MATERIAL 

USED  IN  SAMPLE  CALCULATIONS 

[ten.  denotes  tension;  comp,  denotes  compression;  and 
superscript  a  indicates  allowably 


Symbol 

Value  in 

SI  units 

Value  in  U.S. 
Customary  units 

Density  and  elastic  properties 

p 

1520  kg/m3 

0.055  Ibm/inS 

El 

145  GPa 

21  X  10®  psi 

E2 

16.5  GPa 

2.39  X  10®  psi 

^12 

4.48  GPa 

0.65  X  10®  psi 

^12 

0.314 

0.314 

M21 

0.037 

0.037 

Allowable  stresses 

C7^  (ten.) 

1240  MPa 

180  ksi 

a|  (comp.) 

-1240  MPa 

-180  ksi 

o-a  (ten.) 

55.2  MPa 

8  ksi 

(comp.) 

-207  MPa 

-30  ksi 

ja 

12 

82.7  MPa 

12  ksi 

Allowable  strains 

ea  (ten.) 

0.0087 

0.0087 

(comp.) 

-0.0087 

-0.0087 

e|  (ten.) 

0.00475 

0.00475 

e|  (comp.) 

-0.01764 

-0.01764 

ya 

^12 

0.01846 

0.01846 

TABLE  V.-  PROPERTIES  OF  ALUMINUM  USED 
IN  SAMPLE  CALCULATIONS 


Material  properties 

Symbol 

Value  in 

SI  units 

Value  in  U.S. 
Customary  units 

p 

2780  kg/m^ 

0.101  Ibm/in^ 

E 

72.4  GPa 

10.5  X  10®  psi 

G 

27.4  GPa 

3.98  X  10®  psi 

M 

0.32 

0.32 

T-a 

221  MPa 

32  ksi 

TABLE  VI.-  LIMITS  ON  DESIGN  VARIABLES 


Dimension 

Hat- stiffened 
panel 

Corrugated 

panel 

cm 

in. 

cm 

in. 

t3(0^),  maximum 

0,318 

0.125 

0.318 

0.125 

14(0°),  maximum 

.635 

.250 

Not  applicable 

bj,  minimum 

2.03 

.80 

2.03 

,80 

b2,  minimum 

1.02 

.40 

1.02 

.40 

b3,  minimum 

2.03 

.80 

2.03 

.80 

b4,  minimum 

_  . 

1.02 

.40 

Not  applicable 

NASA- Langley,  1976  L-10770 
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